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We prove a theorem, using the density functional approach and relying on a classical result by 
Lieb and Simon on Thomas-Fermi model, showing that in the thermodynamic limit bulk matter is 
at most semiclassical and coherence preserving. The connection between quantum fluid dynamics 
and density functional theory in the formulation due to Kohn and Sham plays a significant role 
leading to a Vlasov-Poisson system of equations for the Wigner function. Coherence stability is 
achieved by noting that small oscillations in bulk matter are damped by Landau damping. In some 
conditions the initial Wigner function could generate an opposite effect and coherence stability can 
be lost involving higher order quantum effects for a macroscopic body. 

PACS numbers: 71.10.-w, 31.15.Bs, 71.90.+q 

I. INTRODUCTION 

The proof of the stability of the matter, initially due to Dyson and Lennard 1 '- with a mathematical tour de force, 
and largely improved by an elegant argument due to Lieb and Thirring 3 , is surely one of the most important results 
about the existence of bulk matter. The situation is far less clear about the coherence properties of matter, the aim 
of this paper. 

The problem of addressing a many-body system from a quantum mechanical standpoint is rather demanding and 
difficult to treat unless the number of components is really small. A smart approach to this problem arose initially 
with the pioneering work by Thomas and Fermi^ where was shown how density is all one needs to address the 
physics of large atoms, then theorems by Kohn and Hohenberg 6 and Kohn and Sham- made this approach mature 
to analyze the behavior of large systems relying just on the computation of the density. Striking agreement between 
experimental and computed data is so achieved by this so called density functional theory (DFT). 

One of the main drawbacks of this approach, notwithstanding the great successes, should be identified with the 
impossibility to analyze the coherence properties of bulk matter. A question like "is bulk matter a classical or a 
quantum object?" or "does bulk matter preserve coherence?" cannot be addressed just considering density. 

Indeed, these questions can be addressed inside the same theory by a set of results recently obtained^. A Wigner 
function can be associated to a density functional and the coherence properties of bulk matter can be analyzed in the 
thermodynamic limit. 

A relevant result about a many-body quantum theory was obtained by Lieb and Simon 9,10 giving to Thomas-Fermi 
theory a role of limit theory for the quantum counterpart when the number of particles is taken to go to infinity, the 
limit that is of interest for our aims. Thomas-Fermi theory also played a relevant role in the argument of Lieb and 
Thirring on stability of matter. Then, this theory should be considered a substantial aspect of many-body physics. 

Being Thomas-Fermi theory a semiclassical theory, indeed is the leading order of a Wigner-Kirkwood expansion 
of the Green function of the quantum systemAi, we can conclude that bulk matter, in the limit of a large number 
of particles, has classical properties that can be preserved during time evolution. This conclusion about stability is 
obtained by analyzing the time evolution of the Wigner function obtained by time dependent DFT plus the Lieb and 
Simon theorem. Coherence stability means that small perturbations on the system are damped. This is a property of 
the Vlasov-Poisson equations that drive the system in this caseAi. On the other hand, if matter is initially properly 
prepared, coherence stability can be lost and small perturbations can be amplified. Then, classicality, a common 
property, can be lost. 

The paper is so structured. In sec[H] we introduce the model and, by using the semiclassical Wigner-Kirkwood 
expansion and the Lieb and Simon theorem, we show that the limit of the quantum theory for the number of particles 
going to infinity and Ti — > coincide with the Thomas-Fermi model. In order to study coherence properties of bulk 
matter we need a time-dependent model. This is obtained in sec lIIII by the Runge and Gross theorem giving the 
time-dependent density functional theory that translate into a hydrodynamic model for the Thomas-Fermi model. 
Having a hydrodynamic model, a kinetic equation for the Wigner function is obtained in sec II VI and, by a recently 
proved theorem, this is equivalent to the classical limit plus a quantum correction giving the full set of equations 
known as Vlasov-Poisson equations. This set can produce for small perturbations a damping granting coherence 
stability. Exotic initial Wigner functions can generate instabilities. In sec[V]the conclusions are given. 



2 



II. THOMAS-FERMI THEORY AS SEMICLASSICAL LIMIT OF A QUANTUM THEORY 



All the low energy phenomenology is fully described by the Hamiltoniax 
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that is able to describe all the properties of matter as far as we know. We have put N e the number of electrons, JVj 
the number of positive ions, e the electron charge, Z a the number of positive charges for each ion. 

Introducing the density and neglecting for our aims the repulsion between positive ions assumed to be pointlike 
and classical, the quantum system we are going to consider is rather well-known, being the starting point of a lot of 
low energy condensed matter analysis^, having the following energy functional: 
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being T the kinetic energy given by 
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being to the electron mass, p(x) the electronic density computed by the Slater determinant if the Hartree-Fock ap- 
proximation is invoked, as J d 2 X2d 3 xs . . . d 3 XAr^/*(x, x 2 , . . . , xpj)$(x, x 2 , . . . , xn), Z is the number of positive charges 
and e the electron charge. We assume neutrality, that is, the number of electrons, N — J p(x)d 3 x, is the same as 
the number of positive charges in the system. One could transform the functional £ into a density functional £[p(x)], 
bypassing self-consistent Hartree-Fock equations for the Hartree-Fock approximation, if it would be possible to obtain 
the kinetic energy as a functional of p(x). 

We would like to accomplish such a task in the classical limit ?i — > 0. We can reach our aim by referring to standard 
results in nuclear physicsii and Bose condensationM as also given inM: The density matrix can be developed into 
a h — > series, derived from the Wigner-Kirkwood series for the Green function for a single particle moving in the 
Hartree-Fock potential V— 
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whose leading order is just the Thomas-Fermi approximation 
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giving the density of kinetic energy as 
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So, we can conclude that the Thomas-Fermi approximation can represent classical objects being the leading order of 
a semiclassical approximation. 

For our aims this is not enough as we are interested in the limit of a large number of particles in the energy 
functional we started with. Indeed, there is another theorem due to Lieb and SimonSiifl that gives an answer to this 
question: The limit of number of particles N going to inGnity for energy functional |2J) is the Thomas-Fermi model. 
This means that in this limit we recover again the leading order of a semiclassical approximation. So, classical objects 
can be obtained when the number of particles in a system increases without bound because the limit of a number of 
particles going to inGnity coincides with the semiclassical limit Ti — » 0. This is the main result of this section. 

One can ask at this point if this property is stable, that is, if classicality is indeed preserved during time evolution 
or if some instabilities may set in moving the system toward some generic quantum state. In the following we will 
answer to this question relying on the results of this section but to do this we need to study the time evolution of the 
quantum system. We develop this analysis in the next section. 
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III. DFT AND QUANTUM FLUID DYNAMICS 



In the preceding section we have seen that the energy functional can be turned into a density functional. Indeed, 
Thomas-Fermi theory turns out to be the simplest of density functional theories and the relevance of this approach 
relies on a couple of theorems proved by Kohn and Hohenberg&. On the basis of these results the density p(x) contains 
all the informations about the ground state and there is no need to have the full wave function to get the properties 
of the system (first theorem) . Once one has a density functional, a variational principle gives the equations for the 
computation of the density (second theorem). This is generally done numerically by solving a set of self-consistent 
single particle equations as devised by Kohn and Shaml. Hohenberg-Kohn theorems have been extended to the time- 
dependent case by Runge and Grossi^ but now the variational principle changes into a stationary point principle for 
the functional 

A(t)= I 1 ' (V(t)\ihd/dt- H(t)\V(t)). (7) 
J to 

All these results apply straightforwardly to the simplest case of the Thomas-Fermi approximation producing a 
hydrodynamic modelii that can be cast into the form 
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To get this equations we have used the ansatz for the signle particle Kohn-Sham orbital cf>i(x,t) = 
0i(x, i) exp \i— %(x, t)] where we have assumed a single phase for each orbital. As we will see this approximation 
simplifies the computations but does not invalidate our conclusions. 

We can now give explicitly the above set of equations through the density and the velocity field. One has 

d\x 1 

m-^+m(u- V)u + W+-VP F = (12) 
at p 

being V the potential solution of the Poisson equation 

A 2 F(x, t) = 47re 2 [p(x, *) - Z6(x)} (13) 



and Pp the Fermi pressure given by 



P F (x,t) = |-(37r 2 )lp(x,t)l. (14) 

OTO 



that contains all the quantum information. We have obtained a time-dependent Thomas-Fermi model expressed 
through Euler hydrodynamic equations. This set of equations can be transformed into a single Schrodinger equation^ 
but we do not pursue this matter here. 

It is important to notice that the simple knowledge of the density it is not enough to get informations on the 
coherence properties of a many-body quantum system. But, having introduced hydrodynamic equations by the 
density, we can associate to it a Wigner function W(x, p, t) in the phase space containing all the informations about 
the evolution of the system. It is important to notice at this point that, because the Thomas-Fermi model corresponds 
to Euler equations, no collisional term does enter into the kinetic equation for the Wigner function. This opens the 
way to the main result of this paper, that is, classicality, when emerges, is a stable property of a many-body quantum 
system being the system, in the limit of a large number of particles, driven by a classical set of kinetic equations, 
known in the current literature as the Vlasov-Poisson equations. 
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IV. QUANTUM FLUID DYNAMICS AND WIGNER FUNCTION 

For our aim, the set of hydrodynamic equations can be cast in the form 

^ + V.(/m)=0 (15) 

— + (u- V)u + -V(V + £ F ) = (16) 
ot m 

A 2 V(x,i) = 47re 2 [p(x,i) - Z6(x)] (17) 

being Bp = ^-(37r 2 )i/?(x, t)i the Fermi energy. By a recently proved theorem due to Jungel and Wangi& this set of 
equations derives from the self-consistent set of Wigner-Poisson equations 
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p(*,t) = j W(x,v,t)d a v (20) 

being W(x, v, t) the Wigner function. This last set of equations is all one needs to describe low energy phenomenology 
of bulk matter in the limit of number of particles going to infinity. The set we have obtained is the Vlasov-Poisson 
set representing the leading order of a semiclassical approximation to the full Wigner equation 18 . We note that we 
have to keep a quantum term due to the Pauli principle, that is, the fermionic nature of the matter. This shows one 
of the main results of the paper, that is, when the number of particles goes to infinity the coherence is ruled by a 
Vlasov-Poisson set of equations representing the classical limit of the Wigner equation but with a quantum correcting 
term. The correcting term is just the Fermi energy that for a macroscopic body is negligible with respect to the 
self-consistent potential. A similar approximation is generally assumed in Bose condensates for the Thomas-Fermi 
approximation to the Gross-Pitaeskvii equation^. 

This permits us to get finally the proof of our theorem: classicality as an emerging property of a macroscopic 
body is a stable property. This can be seen easily assuming an initial Wigner function, as the one derived from the 
Fermi-Dirac distribution at zero temperature, that warrants Landau damping of small perturbations^ and taking 
into account neutrality of the matter. It is interesting to note that some possible initial conditions may generate 
unstable macroscopic bodies having quantum properties as damping can change into an exponential increase of the 
perturbation and stability could just be reached by properly including higher order quantum effects. 

In order to see the above conclusions at work, we analyze the case of an isotropic superconductor. This implies a 
constant wavefunction, when no current flows, that in turn means that the initial Wigner function will be proportional 
to S(v), being v the modulus of the velocity. One can prove from the theory of Landau damping^ that the condition 
for a small initial oscillation to die out is 
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being Wo(k,v) the initial Wigner function after a Fourier transform in space. So in this case the Landau damping 
factor is null. Indeed, this initial Wigner function is well-known in literature for the Vlasov-Poisson equations giving 
rise to instability 20 and so, in this, higher order quantum corrections are needed to make the system stable. On the 
other side, if we take a Fermi distribution function in the limit of high temperature, this reduces to a Boltzmann 
distribution giving rise to Landau damping proper to a stable classical system. 



V. CONCLUSIONS 



By a sequence of known results we have shown how coherence can be considered for a many-body quantum system. 
Classicality proves to be an emerging property in the limit of a large number of particles as this limit is given by the 
Thomas-Fermi model as also happens to the semiclassical limit h — > at the leading order. Evolution in time of the 
density functional, provided by Runge and Gross theorem translates into a set of hydrodynamic equations for the 
Thomas-Fermi model. These equations are equivalent to the set of Vlasov-Poisson kinetic equations for the Wigner 
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function, corresponding to the classical limit plus a quantum correction given by the Fermi energy, that gives rise to 
Landau damping of small perturbations assuring stability of classicality in the considered limits. 
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